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Abstract 


Assigning available resources to fire stations is a main task of fire depart- 
ment’s administrator in a city. The importance of this problem increases 
when the number of available resources are inadequate. In this situation, 
the goal is to assign the limited available resources to fire stations such that 
the associated penalties of the shortages are minimized. Here, we first give 
a mathematical approach to consider some penalties for the shortage. Next, 
we give an integer program to minimize the sum of associated penalties. The 
proposed model can be used in many other problems arisen from health ser- 
vices, emergency management, and so on. We also propose a heuristic to 
efficiently solve the problem in a reasonable time. Our proposed heuristic 
has two phases. In the first phase, using a greedy approach, our proposed 
heuristic constructs a proper feasible solution. Next, in the second phase, we 
propose a local search to improve the quality of the solution constructed in 
the first phase. To show the efficiency of our proposed heuristic, we com- 
pare our proposed heuristic with CPLEX based on the running time and the 
quality of obtained solutions on two groups of problems (real-word problems 
and randomly generated problems). The numerical results show that on the 
80% of benchmark problems, the obtained solution is the same as CPLEX’s 
solutions. Also, the running time of our proposed algorithm is almost 10 
times better than CPLEX’s running time, in average. 


Keywords: Resource assignment problem; Fire stations; Shortage; Integer 
programming; Heuristics. 


*Corresponding author 

Received 22 January 2017; revised 29 March 2017; accepted 23 April 2017 

E.S. Alavi 

Department of Applied Mathematics, Faculty of Mathematical science, Ferdowsi University 
of Mashhad, Mashhad, Iran. e-mail: effatsadat.alavi@gmail.com 


R. Ghanbari 
Department of Applied Mathematics, Faculty of Mathematical science, Ferdowsi University 
of Mashhad, Mashhad, Iran. e-mail: rghanbari@Qum.ac.ir 


129 


130 E.S. Alavi and R. Ghanbari 


1 Introduction 


There are many decision making problems for a fire department’s administra- 
tor or even for a fire district’s administrator (see [5,9,{C3]; also, (1,3,{7,101,12]). 
One of these decision problems is to assign available resources (e.g., equip- 
ments or fire engines etc.) to fire stations (see [[2], page 524, to find some 
references). In the real world situations, it is possible that there are not 
sufficient resources to satisfy all the requirements of each fire station. Thus, 
some of fire stations may lack some resources to extinguish or to service 
the fire accidents or events in their regions (other nearly problems in emer- 
gency management are studied by researchers, e.g., [2,6,8]). To decrease 
the influences of the shortage for a type of resource, some penalties are con- 
sidered, based on its importance in each fire station. In a fire station, to 
compute the importance of a required resource, we estimate the possibility 
of the occurring of each event in the fire station’s region. Next, based on 
the predefined impact of the events and the estimated possibility for arising 
of the events, the corresponding penalties are considered. In addition, when 
fire incidents occur or rescue services are needed, the responsible fire station 
usually requests additional resources or even firefighters from its neighboring 
fire stations (neighboring services). However, the aim of the problem is to 
assign the available resources to fire stations such that the sum of the corre- 
sponding penalties is minimized. 

Here, we give an integer linear formulation to model the problem. We 
also propose a heuristic to efficiently solve the problem. We finally compare 
the performance of our algorithms with ILOG OPL(CPLEX solver)’s result 
on some randomly generated test problem and some real world problems. 

The remainder of the paper is organized as follows: In section 2, we give 
an integer program to model the problem. The proposed algorithm is pre- 
sented in section Bl In section B, we show the efficiency of our proposed 
algorithm in comparison with ILOG OPL on some randomly generated test 
problems and some real world problems. We conclude in section BH. 


2 Mathematical model 


We here consider all resources are independent. It is obvious that when two or 
more types of resources are dependent, one can consider them as a package. 
Therefore, we can assume the shortage exists for a resource. It is obvious 
that the model and related solution approach can be used as a procedure 
within a successive approach for solving the problem with shortage in two 
or more types of resources. As mentioned in section [], before providing the 
integer program, we must give some necessary estimation or computation 
to determine the corresponding penalty of the recourse’s shortage in a fire 
station. Table [I] presents some parameters or given data of the problem. 
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Table 1: Parameters or given data 


the availability of the resource. 


gq: the number of events. 

n: the number of fire stations. 

tr: the required number of the resource to service the event k. 

Dkj: the frequency of the event & in the region of station 7. 

prrj: the possibility for occurring the event k in the region of station 7. 
Wr! the importance (weight) of the event k. 

N(j): the set of the fire station j’s neighboring fire stations. 

By: the importance (weight) of the fire station j. 


Note 1. In Table, pr.,., w., and £. are in [0, 1]. 


Note that pz; is obtained previously. In the other words, we can consider 
the frequency of the event k in region of the fire station 7 for example for 
a period of ten years. But, pr; is given by managers, and it shows the 
possibility of occurring the event k in the region of fire station 7. To more 
clarify, we give an example. Consider, a school is recently founded in the 
region of the fire station 7, and suppose there was not any school in this region. 
Therefore, some events have rarely occurred there. But, the possibility of rare 
occurring or never occurring events is now raised. Thus pr;; and pp; can be 
completely different. When all the situations are the same, 3; imposes the 
manager’s decision to assign the resource to the fire station 7. As an example, 
consider there are two fire stations with the same conditions, need a type of 
equipment. Also, consider one of them is in the center of the city, and the 
other is on the outskirts. Suppose that there is one equipment. In this 
situation, it is usually preferred to assign this equipment to the fire station 
that is in the center of the city. Therefore, the manager can consider the 
higher weight to the fire station that is in the center of the city to impose 
her/his preferences in the model. 


To compute the associated penalty, we must do some necessary computa- 
tions. We consider the relative frequency for the event / in the region of the 
fire station j as follows: 


Ukj = _ , (1) 
Do Pj 
j=l 


Remark 1. Note that one can consider u,; introduced in (MM) as an esti- 
mation of the probability of occurring the event k in the region of the fire 
station 7. 


Definition 1. Let 
Ukj = max{ux;, PTR; }, (2) 
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where pr; is given by the manager and ux; is introduced in (fl). Now, we 
define vz ;w, as the impact of the event k on the station j, for 7 = 1,...,n 
and k =1,...,q. 

Note that we could not find any definition for vz; in the literature since 
this model was proposed for a real problem that occurred in Mashhad. We 
did not have no more information in this real-world situation. Therefore, 
using experts in Fire Department of Mashhad, we proposed this formula. In 
this formula, we consider historical data pz; and experts consideration pr,;. 
Thus, w,v,%; shows how important the event k in the region of the station j. 


Definition 2. We denote the possible events in the station j by Q,; de- 
fined as follows: 


Qi = {kl fp > 0,03 > Of. (3) 
where v;; is defined in (2) and f;, is defined in Table I. 
We now introduce the decision variables in Table 2. 


Table 2: Decision variables 


x;: shows the number of the resources that are assigned to station 7. 

Yjk: Shows the shortage in the station 7 to service the event k. 

hjx: shows the shortage in the station 7 to service the event k with considering 
the neighboring services. 

dj,: if station j can not service the event k, takes 1, otherwise takes 0. 


2.1 The objective function 


The objective of the problem is to minimize the sum of the associated penal- 
ties. We here consider three penalties for servicing the event k in the station 
j based on the lack of ability for providing responsible service, the number of 
shortages and the number of shortages with considering neighboring services. 
Therefore, we consider the following penalty for the station j: 


> (djx + Yjk + hjk) WkVE; (4) 
kEQ; 


where Q,; and vz; are defined in (@) and (Q), respectively, djx, yjx, and hjx 
are defined in Table 2, and wz, is defined in Table [1 To explain (@), we give 
two illustrative examples. 


Illustrative Example 1. Consider there are five fire stations and only one 
resource exists. Also, the structure of the neighboring fire stations is given 
in Figure J. Suppose that each of the stations 1, 3, 4, and 5 need one unit 
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C1 2 me 3% 


Figure 1: the graph of the neighborhoods in the illustrative Example 1. 


of the resource and all situation are in the same. From Figure [I], the station 
3 is the best candidate, because, the station 3 while servicing the occurred 
events in its own region, can participate to service the occurred events in the 
stations 4 and 5’s regions. One can easily obtain the same result by using 
the sum of the penalties defined in (A) for each fire station. 


Illustrative Example 2. Consider there are two fire stations and three 
possible events. Also, each fire station is without a neighbor (i.e. N(1) = 
N(2) = 9). The other data are given as follows (see tables IH2 and equations 
(2) and @)): 


T=8, Qi = {2,3}, Qe = {1,2,3}, 
fe (0.2, 0.3, 0.7], 


| + 


[4, 6, 2], w= 


oO 


7 
» — [0-2 0.9 0.7 
~ |0.4 0.9 0.2 


Now, we give the motivation for considering d;,; in our proposed model. If 
d,; is removed from (@), we have the following penalty for the station j: 


S> (yin + hj) WeeR;- (5) 
keEQ; 


Since N(1) = N(2) = 9, (8) can be considered as follows: 


S> (vik) War R;- (6) 


kEQ; 


Now, consider two feasible solutions, (a): x1 = 2% = 4, (b): x, = 2 and 
2 = 6. For (a), we have 


yi2 = 2, yi3 = 0, yor = 0, Yoo = 2, yo3 = 0, 


2x0.9x0.38+2 x 0.9 x 0.3 
and based on (@), the sum of penalties is . — a = 


1.08 0.7 
a7: For (b), we have 


yi2 = 4, yi3 =0, yor = 9, Yoo = 9, yo3 = 0, 
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4x0.9x 0.3 1.08 


and based on (@), the sum of penalties is . Thus, when 


we do not consider d,; in the model, there is no difference between (a) and 
(b). But, in (b), only the event 2 in the station 1 is not covered. Whereas 
in (a), the event 2 in stations 1 and 2 is not covered. Now, if we consider (Al) 
as the associated penalty, then (b) is preferred to (a). 


2.2 The model 


Now, the proposed integer program is given as follows. 


min z= > B; S- (djn + Yjk + hye) WeVE; (7) 
j=l kEQ; 
s.t. 
j=l 
Viele Nie cg PODS (9) 
hjr = fe —@j — S- U5 Vj €{1,...,n} ke Q;. (10) 
JEN(3) 
djk < Yjr Vj €{1,...,n} ke Q;. (11) 


dix > 5, Vj €{1,...,n} ke Qj, where Mp = max fy. (12) 
f k=1,...5¢ 


Ujk» hy € NU {0}, and dix e {0,1} VIE {1,...,n} ke Qj. 
x; © NU {0} Vj € {1,...,n}. 


The constraint (8) shows that the number of assigned resources must be less 
than or equal to 7. From Table 2, yz; shows the shortage for servicing the 
event k in the station 7. Thus, if we assign 2; unit of the resource to station 
j, then yx; is max{0, f, — x;}. Based on constraints (Q), if «; > f;,, then 
the corresponding constraint is redundant. Since the objective function is 
minimization, then yz; takes zero. In the other case, if x; < fx, then the 
corresponding constraint ensures that y; > f, — x; > 0. Similar to above, 
since the objective function is minimization, then yz; takes f,—x,;. Therefore, 
constraints (@) ensure that y;,; is max{0, f;,—x;}. Similar to constraints (9), 
constraints (TM) ensure that hj = max{0, fr — 23 — jen (j) tj}. Together 
constraints (I) and (12) guarantee, if y,; is zero, then d;; is zero and if 
Yrj > 0, then dz; is 1. 
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Algorithm 1 Greedy phase 


Step 1 {Initialization} Let x* = (0,...,0)2,,, and let z* = z(a*). 


Step 2 {Main loop} 
For t= 1 to T do: 
- Let z, = +00. 
For j = 1 to n do: 
- Lett @=2*, & =2@;+1 and z= 2(Z). 
- If Z < z,, then let z, = Z and 7 = j. 
endfor 
- Let x; = vs +1, and let z* = z,. If z* = 0, then stop, the optimal 
solution is x*. 


endfor 


Step 3 Return 2” as an initial solution. 


3 A heuristic algorithm 


In this section, we propose a heuristic algorithm to solve the problem. Our 
proposed algorithm has two phases. In the first phase, we propose a greedy 
procedure to construct a feasible solution (the greedy phase). Then, in the 
second phase, we apply a neighborhood search on the constructed solution 
(the improvement phase). 


3.1 Greedy phase 


In this phase, we construct a feasible solution. At first, we let 71 =--- = 
Ly = 0. Then, in each iteration, we compute the corresponding difference of 
the objective function for each fire station, when the number of the assigned 
resources in a fire station is increased one unit. Then, we select the station 
j to increase x;, when the corresponding difference of the objective function 
is the best. 

To evaluate x = (71,...,%n)7, let yz; = max{0, fy — xj}, and let hj, = 
max{0, fi — 27 — EN (s) x7}. Also, if yx; is zero, then let dy; = 0 and 
otherwise let dy; = 1. Now, the corresponding objective value of x can be 
computed by (1). We denote the corresponding objective function value of x 
by z(xz). The greedy phase is given in Algorithm [1]. 
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Note 2. In Algorithm [, if we can find a solution 2«* with z(a*) = 0, then 
we have found a solution without any penalty. Therefore, we terminate the 
algorithm, and we do not use the phase 2. 


3.2 Improvement phase 


Here, we describe the neighborhood search procedure. We use a simple local 
search to improve the quality of the solution constructed by Algorithm []. We 
first give some definitions. 


Definition 3. We denote the maximum number of requirements in the 
station 7 by b; and define it as follows: 


Bp ely i a yhsem (13) 


where Q, is defined by (QB). 


It is evident that if x; > b;, then we do not have any shortage in the 
station 7. 


Definition 4. Let 


a min{ f, — aj|k © Q;, fr > a;}, v5 < b,, 
A; _ 6 v5 = b;, (14) 


where Q, is defined by (@). A; is the minimum number of the resources that 
if assigned to station j, an uncovered event (event that does not receive a 
response service) will be covered. 


To explain our proposed local search, we need to define a move in the 
search space (see [I]). For the station j with x; < b;, we define 


M; = {jlaz > Aj. 5 #3}; (15) 


where A, is defined by (14). Note that, M; is the set of stations that their 
current availability is more than the minimum requirement of the station 7. 
Now, for the station 7 with x; < bj, a feasible solution can be obtained from 
the current feasible solution by x; = 7;+Aj; and xj = xj—A,j, where j € M;. 
We denote this move by Mov(j, 7, A;). 


Note 3. When a move is done, an effective procedure to compute the dif- 
ference of the objective function value, can decrease the computational time. 
Consider x is a feasible solution and % is the obtained feasible solution by 
Mov(j,j, Aj). Also, consider z(z) and z(Z) are the objective function values 
of « and Z, respectively. We now define 
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Algorithm 2 Local search (the improvement phase). 


Step 1 Let 2* is the solution constructed by Algorithm I], and let 2* = z(x*). 
Step 2 Let S = {j|xj < bj}. 
Step 3 For j = 1,...,|S| do: 


3-1 Compute A; by (14). 
3-2 For j € Mj, defined by (Gh), if A¥; , <0, defined by (4G), then 
Mov(j,j, Aj), update x* and z*, and goto Step 2. 


Step 4 Return <* as an approximate solution. 


AG 5.A; = 2(Z) — 2(x), (16) 


where 2(.) is defined by (@). It is obvious that AZ; , can be computed by 


subtracting the sum of increased penalties of station j and its neighbors from 
the sum of the decreased penalties of station 7 and its neighboring stations. 
Therefore, the other stations don not play any role in the computation of 

Vi Aj} thus in the implementation, the value of At 5, A, can be computed 


effectively. Also, One can use an approximation of A; = 


in the implemen- 
3,4; e impleme 


tation (similar to [A]). 


We now give our proposed local search (as the improvement phase) in 
Algorithm 2. 


Proposition 1. Our proposed heuristic gives a feasible solution. 


Proof. In the greedy phase (Algorithm I), our proposed heuristic begins from 
zero (a = 0), and then in each iteration of Step 2, only one unit is added to 
one entry of x. Thus, the return solution is a non-negative integer vector. 
The second phase starts by the feasible solution given by phase 1. In each 
iteration of Step 3 in Algorithm 2, each move Mov(j,j,Aj;) is done when 
x; => A; (see (14)) and A; > 0 defined in (14) is an integer number. Thus, 
after each move the solution remains feasible. 


Note 4. There are some problems such that our proposed heuristic does not 
give the optimal solution (e.g., MASH2 in table B). Also, we are not able 
to prove that our proposed heuristic is an ¢-approximation algorithm [1]. 
However, the numerical results show that our proposed heuristic gives very 
good solution in an acceptable time for real-world problems (see Table B). 
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4 Numerical experiments 


To show the efficiency of our proposed algorithm, we compared the run- 
ning time and the quality of the arrived solution of the proposed algorithm 
with ILOG OPL 6.3 (CPLEX 12.1) on two groups of test problems. In the 
first group, we used problems given by Fire Department of Masshad (MASH 
problems in Table B)). In the second group, we generated some test problems 
randomly (RAND problems in Table B)". We implemented our algorithm in 
Visual Studio C++ 2010. All implementations were done on a Notebook i5 
2.5GHz with 4GB of RAM when only one core was used for our proposed 
algorithm and all cores were used for OPL. All results are given in Table B. 


Note 5. In some test problems, OPL terminated before convergence be- 
cause the memory was not enough. In these situations, we reported the 
truncated results of OPL. 


In Table B, Gap is computed as follows: 


Z-2 
Gap = ——, 
z 


where z is the objective value obtained from our algorithm and 7 is either the 
objective value of the best integer solution found by OPL or the objective 
value of the best node reported by OPL. 


The mean of gaps of our proposed algorithm and the best integer solution 
found by OPL is only 3.7e-4._ Despite OPL used all cores to perform its 
computations and the implementation of our algorithm used only one core, 
the running time of our proposed algorithm is very less than OPL’s running 
time. From table B, it is clear that, in 63% of problems, the given solution 
is equal to the best integer found by CPLEX. In 17% of problems, the given 
solution is better than CPLEX’s solution. In these problems, CPLEX was 
terminated because the memory was not enough. In 20% of problems, the 
given solution is worse than CPLEX’s solution. But in these problems the 
maximum gap from feasible solution given by CPLEX is only 9.6e—3 and the 
running time of our proposed algorithm was 20 times better than CPLEX’s 
running time, approximately. In some problems, CPLEX is not able to find 
the optimal solution. In these situations, the lower bound given by CPLEX 
is suitable for comparing the result. From Table 8, the maximum gap from 
the lower bound is 1.2e — 2. Thus we can say that the maximum error of our 
heuristic on problems is less than 2e — 2 when the running time is almost 10 
times better than CPLEX’s running time. Therefore, from Table Bl, we can 
conclude that our proposed algorithm is applicable for solving the medium 
scale problems in a real time. 


1 For acquiring hold of the test problems, please contact the corresponding author (rghan- 
bari@um.ac.ir) or download at http://profsite.um.ac.ir/~rghanbari/IJNAO.rar 
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Table 3: Numerical Results 
Our Algorithm OPL Gap 

Problem n qd Best Time (sec) __Best(integer) _ Best(node) _Time(sec) __ Integer Node 
MASH1 34 20 883.5649 0.16 883.5649 883.5649 65.52 0 0 
MASH2 34 20 = 3100.1002 1.03 3098.9425 3074.4241 1729.29 de-4 8.4e-4 
MASH3 34 20 731.2355 0.67 731.2355 731.2355 267.40 0 0 
MASHA4 34 20 6538.2014 2.87 6536.7751 6509.6999 1012.93 2e-4 4.4e-4 
MASH5 34 20 591.8589 0.83 591.8589 591.8589 141.85 0 0 
MASH6 34 20 630.4577 0.05 614.9958 614.9958 5.10 9.6e-3 9.6e-3 
MASH7 34 20 0 0.00 0 0 0.08 0 0 
MASH8 34 20 2164.0976 1.28 2164.0976 2147.2527 1285.28 0 7.8e-4 
MASH9 34 20 530.2508 0.41 530.2508 530.2508 53 0 0 
MASH10 34 20 841.3675 2.87 841.3675 791.8676 873.59 0 2.76¢e-2 
MASH11 34 20 222.7586 1.20 222.7586 222.7586 26.21 0 0 
MASH12 34 20 7034.8948 3.62 7034.8948 6976.9721 689.24 0 8.3e-3 
MASH13 34 20 929.6792 0.27 929.6792 929.6792 23 0 0 
MASH14 34 20 815.4193 1.57 816.1045 756.3611 1189.90 -4e-4 3.36e-2 
RAND1 5 00 89.1726 0.00 89.1726 89.1726 0.03 0 0 
RAND2 30 00 4766.9702 0.08 4766.9702 4766.9702 89 0 0 
RAND3 60 00 9055.7206 11.79 9055.7206 9042.4482 698.93 0 1.5e-3 
RAND4 100 00 = 7181.3464 87.25 7181.5357 7123.6819 895.88 -3e-4 8.1e-3 
RAND5 5 50 = 2492.4775 0.00 2492.4775 2492.4775 4.35 0 0 
RAND6 30 50 -5116.5434 1.28 5118.1958 5075.0990 850.25 -le-4 8.2e-3 
RAND7 60 50  3993.0682 16.35 3993.0682 3948.5215 784.03 0 1.1e-2 
RAND8 100 50 —9308.8906 66.46 9309.9365 9047.4370 45.04 -le-4 2.9e-2 
RANDY 5 200 168.5741 0.00 168.5741 168.5741 0.11 0 0 
RAND10 30 200 = 933.2573 0.66 933.2573 933.2573 154.83 0 0 
RAND11 60 200 3437.8746 15.98 3439.9514 3389.9545 706.86 -6e-4 1.4e-2 
RAND12 100 200 3186.5378 131.57 3185.4674 3159.9767 762.33 3e-4 8.4e-3 
RAND13 5 250 22.2134 0.00 22.2134 22.2134 0.2 0 0 
RAND14 30 250 93.9136 1.36 93.9136 93.9136 9.44 0 0 
RAND15 60 250 1076.3503 12.33 1074.8568 1063.2535 890.00 14e-4 1.2e-2 
RAND16 100 250 2707.3912 123.80 2706.2883 2667.7602 860.47 de-4 1.5e-2 

‘: truncated (Out of memory). 

*: optimal. 


5 Conclusions and further works 


Assigning available resources to fire stations is a main task of fire depart- 
ment’s administrator in a city. The importance of this problem increases 
when the number of available resources are inadequate. In this situation, the 
goal is to assign the limited available resources to fire stations such that the 
associated penalties of the shortages are minimized. Here, we first gave a 
mathematical approach to consider some penalties for the shortage. Next, 
we gave an integer program to minimize the sum of associated penalties. The 
proposed model can be used in many other problems arisen from health ser- 
vices, emergency management, and so on. We also proposed a heuristic to 
efficiently solve the problem in a reasonable time. Our proposed heuristic 
had two phases. In the first phase, using a greedy approach, our proposed 
heuristic constructed a proper feasible solution. Next, in the second phase, 
we proposed a local search to improve the quality of the solution constructed 
in the first phase. To show the efficiency of our proposed heuristic, we com- 
pared our proposed heuristic with CPLEX based on the running time and 
the quality of obtained solutions on two groups of problems (real-word prob- 
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lems and randomly generated problems). The numerical results showed that 
on the 80% of benchmark problems, the obtained solution was the same as 
CPLEX’s solutions. Also, the running time of our proposed algorithm was 
almost 10 times better than CPLEX’s running time, in average. 

Finally, areas of further research within the framework of our proposed model 
and heuristic include 


1. The second phase of our heuristic was a local search; so we expect 
that if we can propose a advanced search method such as tabu search 
or variable neighborhood search and so on, then the quality of the 
obtained solution will be increased. 


2. In the proposed model, we assumed that some parameters were given 
in a numerical form (e.g., py; and prz;). We guess that the reliability 
of model is increased when the value of these parameters are given 
by intervals or even by linguistic variable. It is evident that in this 
situation, some advanced robust optimization method must be used to 
solve the model. 


3. The proposed model is suitable when we have the shortage in our exist- 
ing resources, and this shortage is very far from the minimum required 
resource to cover an event. When the level of shortage is not high or 
there is not a shortage in each resource, we can give a model to cover 
the events when they can be occurred in parallel. 
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